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Abstract 

In the present paper we discuss the modifications introduced into the 
first-order Fermi shock acceleration process due to a finite extent of dif- 
fusive regions near the shock or due to boundary conditions leading to 
an increased particle escape upstream and/or downstream the shock. In 
the considered simple example of the planar shock wave we idealize the 
escape phenomenon by imposing a particle escape boundary at some dis- 
tance from the shock. Presence of such a boundary (or boundaries) leads 
to coupled steepening of the accelerated particle spectrum and decreasing 
of the acceleration time scale. It allows for a semi-quantitative evaluation 
and, in some specific cases, also for modelling of the observed steep par- 
ticle spectra as a result of the first-order Fermi shock acceleration. We 
also note that the particles close to the upper energy cut-off are younger 
than the estimate based on the respective acceleration time scale. In Ap- 
pendix A we present a new time-dependent solution for infinite diffusive 
regions near the shock allowing for different constant diffusion coefficients 
upstream and downstream the shock. 

Keywords: cosmic rays - Fermi acceleration - acceleration time scale - shock 
waves 

1 Introduction 

In the test particle approximation, the first-order Fermi shock acceleration with 
infinitely extended diffusive regions upstream and downstream the shock leads 
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to a power-law particle spectrum with the spectral index 
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and the acceleration time scale 
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where the index "1" ("2") indicates respectively the upstream (downstream) 
quantity, U is the shock velocity, R = U1/U2, « is the spatial diffusion coefficient. 
For a review of the results referring to the diffusive acceleration mechanism one 
should consult some of the numerous review papers (e.g. Drury 1983; Blandford 
& Eichler 1987; Berezhko et al. 1988; Jones & Ellison 1991). 

In discussions of the acceleration process, the background conditions are 
usually considered with particle diffusion coefficients changing at most moder- 
ately with the distance from the shock. As a result, there are infinite diffusive 
regions for cosmic ray particles. Then, particles are removed from the accel- 
eration region near the shock only due to advection with the general plasma 
flow far downstream. However, if the waves responsible for particle scattering 
are created due to the process of cosmic ray streaming instability upstream the 
shock, e.g., the finite amplitude waves resonant with particles near the spec- 
trum cut-off energy could be created only close to the shock within a finite time 
available for the acceleration process. The analogous situation is the case if 
the shock propagates through the finite volume of turbulent plasma. Then, the 
energetic particles diffusing far from the shock will encounter the conditions en- 
abling them to permanently escape from the shock. A finite extent of the shock 
wave to the sides and some particular boundary conditions may also allow for 
such escape. The situation can be qualitatively modelled by introducing the up- 
stream (Berezhko et al. 1988) and/or downstream boundary for the energetic 
particle escape. 

The conditions mentioned above are sometimes discovered in analysing in 
situ measurements near heliospheric shock waves. Then, an attempt to analyse 
the data within the standard approach, involving equations (1.1-2) to describe 
the first-order Fermi acceleration, may fail. For example, basing on such an 
analysis Bialk & Droge (1993) rejected the possibility of the first-order acceler- 
ation at the considered shock wave and suggested the second-order acceleration 
downstream the shock to play a role, instead. Unfortunately, to date there 
is no theory available describing effects of the enhanced particle escape at the 
acceleration process and discussions of the 'difficult' cases have to be based on 
numerical methods. As a step forward, in the present paper we develop a simple 
analytic theory allowing for evaluation of measurements in terms of the diffusive 
length scales for energetic particles. It can become a starting point for more 
elaborate computations or numerical modelling. The theory describes modifi- 
cations introduced into the acceleration process due to finite extent of diffusive 
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Figure 1: The configuration of the shock with the upstream free escape 
boundary ("u.f.e.b.") and the downstream free escape boundary ("d.f.e.b."). 
The situation is presented as seen in the shock rest frame. Between the shock 
and a given escape boundary there is a diffusive layer with some finite diffusion 
coefficient k. Behind the boundary one assumes infinite diffusion coefficient al- 
lowing for a free streaming to infinity (escape) of particles hitting the boundary. 

regions near the shock (Figure 1). Below, in Section 2 we discuss a simple 
time-dependent solution of the diffusion equation for cosmic ray particles at the 
shock with the upstream and/or the downstream escape boundary. To obtain 
analytic solutions, we consider a simple case with constant diffusion coefficients, 
leading in the stationary situation to power-law cosmic ray spectra. It enables 
us to discuss the relation between the spectrum inclination and other param- 
eters of the acceleration process. In Section 3, we derive the particle spectral 
index and the acceleration time scale as a function of the boundary distance 
from the shock. The introduced particle sinks at escape boundaries lead to 
steepening of the spectrum accompanied by a substantial - a factor of two or 
three for reasonable spectral indices - decrease of the acceleration time scale. 
We discuss the dependence of these quantities on the upstream and downstream 
boundary distance from the shock. Finally, in Section 4, we shortly summarize 
the results. The presented theory allows one to interpret the above mentioned 
Bialk & Drogc (1993) data in terms of the first-order acceleration. We also note 
that the particles close to the upper spectrum energy cut-off are younger than 
the estimate based on the respective acceleration time scale valid at infinite 
times. The new time-dependent solution for the particle distribution function 
in the case of infinite boundary distance is presented in Appendix A. It is a 
generalization of the Toptygin (1980) solution for the case of different constant 
diffusion coefficients upstream and downstream the shock. 
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2 On solution of the time-dependent problem 



The problem of time-dependent solution of the diffusion equation for cosmic ray 
particles accelerated at the shock has been discussed by Drury (1983, 1991). 
Below, we will follow the approach described in the former of the mentioned 
papers. In order to produce power-law particle momentum spectra with the 
finite extent of free escape boundaries one must restrict considerations to the 
momentum- independent diffusion coefficient. As in such conditions the spatial 
dependence of k does not lead to any qualitative changes of the acceleration 
process, we consider the most simple situation with the diffusion coefficient spa- 
tially constant, separately upstream and downstream. As mentioned previously, 
we consider the particle acceleration process in the test-particle approximation 
at the non-modified shock structure, i.e. with flow velocities constant outside 
the shock. 

Let us consider a planar shock wave propagating along the rr-axis of the 
reference frame, with the velocity in the upstream plasma frame directed toward 
decreasing x values. In the text below we will refer to the spatial co-ordinates 
in the shock rest frame, where the co-ordinate x upstream the shock is negative. 
We are going to solve the diffusion equation for the isotropic part of the cosmic 
ray distribution / = f(t,x,p): 

where i = 1 or 2 and the mono-energetic source function is taken as 

Q(t,x,p) = Q S(x)5(p-p )H(t) , (2.2) 

where Qo is a constant, 5 is the Dirac delta function and H is the Heaviside 
step function. We look for the solutions satisfying the upstream (x = —L\) and 
downstream (x = L 2 ) free escape boundary conditions 

f(t,-L u p)=0 = f(t,L 2 ,p) , (2.3) 
and the matching conditions at the shock (x = 0): 



[/] = , (2.4) 



dx 3 dp 



-Q 6(p-p )H(t) , (2.5) 



where the square brackets denote the differences between the upstream and 
downstream values of the given quantity across the shock. Following a standard 
approach described by Drury (1983) we start with the Laplace transform with 
respect to time of the distribution /, 
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g(s,x,p)= / e- st f(t,x,p)dt , (2.6) 
Jo 

leading to the following form of Equation (2.1) at p > p : 

The solutions to this linear equation have the form g(x) = C + exp ((3 iy+ x) + 
C_ exp (J3i-x) (i — 1, 2; C± — const). With the definition Tj = 'ini/Uf, the 
exponents /3, are given as 



Pi,± = — 77 (l "F VI + Ti«) 



(2.8) 



Boundary conditions (2.3) for the function g(s,x,p) are g(s,— L\,p) = = 
g(s,L2,p) and the upstream and downstream solutions of Equation (2.7) can 
be written as 



gx(s,x,p) = Ci(s,p) 



e 01 -~ x -e 
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g 2 (s,x,p) = C 2 (s,p) 
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One may note that, contrary to the solutions with no escape boundary, the 
present ones make use of both the indices /?+ and /?_ . The matching conditions 
at the shock for the function g(s,x,p) are derived from Equations (2.4,5) as 
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dx 3 dp 
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By imposing the condition (2.11) at the solutions (2.9,10) we obtain 

" fjln — • Ll e / 3 i, + :c 



ffi(s,a;,p) = g (s,p)- 



(2.11) 
(2.12) 

(2.13) 
(2.14) 
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where go(s,p) = g(s,0,p) is the Laplace transform of the distribution function 
at the shock. Its functional form can be derived with the use of the condition 
(2.12) and Equation (2.8) as 
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g 2 (s,x,p) = g (s,p)- 
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g (s,p)= r J Q ° , ( P A a(S) H(p- P ,) , (2.15) 



where 
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The distribution f(t, x,p) can be formally derived by inverting the respective 
upstream and downstream transforms (2.13) or (2.14): 

f j (t,x,p) = — e st g,(s 1 x 1 p)ds (.7 = 1,2) , (2.17) 

2m J-ioo 

where the path of integration lies to the right of all singularities of the inte- 
grand. The result of such inversion for a particular case of L\ = L 2 = oo, and 
for constant diffusion coefficients which can be chosen independently upstream 
and downstream and not necessarily be equal, K\ ^= n 2 , we present in Appendix 
A. However, even without actually performing the integration one is able to 
extract from Equations (2.13-15) quite a lot of information (Drury 1983). The 
asymptotic behaviour at large times is obtained by looking just at the contribu- 
tion of the rightmost singularity of the integrand, here a simple pole at s = 0. 
It gives the power-law steady state spectrum at the shock, /(oo,0,p) = fo(p), 
with the spectral index a = a(0): 

1 1 e~^~ I 
—+n ZWf • ( 2 - 18 ) 

The first term in the product at the right-hand side represents the spectral index 
a for the shock with infinite diffusive regions (cf. Equation 1.1), and the two 
parts of the second term describe modifications due to particle escape through 
the upstream and downstream boundaries. 

With the use of Equations (2.15,16,18), in the solution (2.17) taken at the 
shock (x = 0) one can separate a part representing the limiting stationary 
solution fo(p) multiplied by a factor describing the time dependence of the full 
solution: 

fo(t,p) = fo(p)- l\{t')dt' , (2.19) 
Jo 

where, with the notation Aa(s) = a(s) — a , 
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Aa(s) 



(2.20) 



One can readily show that the integral 



/ ip(t)dt 



= 1 



(2.21) 
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If we denote 



h(s) = exp [Aa(s) In (p/p )} 
the mean acceleration time from p to p is given as 



(2.22) 



T acc (p, Po )= f tiP(t)dt=^{s = 0) 



(2.23) 



With Equations (2.16, 18.22) and assuming \n(p/p ) — 1 the above formula 
yields the acceleration time scale 



For Li — > oo (i = 1, 2) the expressions (2.18,24) reduce to the standard formulae 



3 The acceleration process with free escape bound- 
aries 

Let us consider variations of the particle acceleration time scale (2.24) and the 
spectral index (2.18) with the escape boundary distance. For the discussion, we 
choose the simplest conditions with K\ = «2 and cither L\ = Li = L, or one 
Li = oo, but the generalization to more general conditions is a straightforward 
one. Below, a boundary distance is expressed in the units of the respective 
diffusive length scale, Ldiff.i = Ki/Ui (i = 1, 2). 

In Figure 2, we present the acceleration time scale variations due to changing 
the escape boundary distance. A substantial decrease of this time scale at small 
L should be noted. The reason for that is clear. The particles which could diffuse 
for a long time in the case with escape boundaries in infinity, for finite L will 
probably disappear from the acceleration process as their chance to wander far 





(2.24) 



(1.1,2). 
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Figure 2: The dependence of the characteristic acceleration time, T acc , normal- 
ized to the acceleration time scale at the shock with infinite diffusive regions, 
T acc ,o, versus the distance L of the free escape boundary (given in units of the 
respective diffusive scale i = 1 or 2). Full lines show the relations for 

the case of equally distant upstream and downstream escape boundary (Li = 
L 2 = L) with different compressions: R — 4 (lower line), R = 3 (intermediate 
line) and R = 2 (upper line). The lines with long dashes represent the above 
relation for downstream particle escape {L\ = oo, = L) and those with short 
dashes for upstream particle escape (Li = L, L 2 — oo). For picture clarity, in 
the last two cases we give curves for our limiting compressions, the lower curve 
for R — 4 and the upper curve for R = 2. 
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Figure 3: The dependence of the stationary particle spectral index a versus 
the distance L of the free escape boundary (given in units of the respective 
diffusive scale Ki/Ui, i=lor 2). The meaning of different curves is the same as 
at Figure 2. Here, the shock compression values are given near the respective 
curves. 



Figure 4: The dependence of the characteristic acceleration time ratio 
Tacc/Taccfi (cf. Figure 2) versus the spectral index (cf. Figure 3). The 
meaning of different curves is the same as at Figure 2, the shock compression 
values are given near the respective curves. 
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off the shock is higher. So, with small L, the spectrum is set up only by particles 
diffusing very close to the shock with a short mean time between the successive 
shock crossings, resulting in the short acceleration time. The dashed lines show 
how the acceleration time is affected by placing the escape boundary upstream 
only or downstream only. Here and in the figures below, one may note an 
apparent asymmetric influence of the upstream and the downstream boundary. 
The difference results from the fact that the mean time between successive 
particle interactions with the shock is related to the diffusive scale different on 
both sides of the shock. In Figure 3, variations of the stationary spectral index 
a® are presented. The expected steepening of the spectrum with decreasing L 
due to the accompanied increased particle escape is presented. For larger L, the 
curves quickly converge to the limiting values 3i?/ (R— 1). Imposing the escape 
boundary downstream only has a much smaller influence on the spectral index 
than the upstream boundary at the same distance (in our units of Ldiff) from 
the shock. This asymmetry can be easily understood by comparing the mean 
times for particles at the shock to reach the escape boundaries. For small Li (i = 
1 or 2), when the particle diffusive streaming dominates over advection, this time 
scale can be estimated as tj w L 2 /fy. For the assumed equal values for upstream 
and downstream diffusion coefficients the respective diffusive distances relate as 
Ldiff, 2 = R Ldiff, l and t 2 = R 2 t\ . So, the escape through the upstream 
boundary will be a factor t%/ti — R 2 greater than the downstream escape, and 
thus it will be mostly responsible for increasing the spectrum inclination. For 
the same reason the downstream boundary will influence the acceleration time 
scale to the higher degree. For larger Li particle advection with the general 
plasma flow modify these estimates, but our conclusions are still true. Always 
the boundary allowing for quick particle escape preferentially acts to increase 
the spectrum inclination, while the other one will provides the lower limit for 
the acceleration time scale. 

A comparison of the spectrum and the acceleration time scale at various L 
is presented at Figure 4. Along the curves, the parameter L decreases to the 
right, i.e. with increasing a. One may note that steepest parts of the curves 
are to the left, near the limiting spectral index. For all the considered shock 
compression values it is possible to decrease T acc by about two times (more for 
larger R, or if the particle escape boundary occurs only downstream) when the 
spectral index steepens on Aa w 1.0 . 

4 Discussion 

The presented results allow one to model a wider range of cosmic ray spectra as 
the output of the first-order Fermi acceleration at shock waves, including dis- 
tributions with a spectral index steeper than the value given by Equation (1.1). 
One should note that accelerating particles with a steep spectrum due to finite 
escape boundaries allows to obtain spectra extended to somewhat higher ener- 
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gies, due to decrease of the acceleration time scale with respect to the one given 
in Equation (1.2). Another point to be mentioned here is an apparent asymme- 
try of the acceleration process with respect to imposing the escape boundary 
upstream or downstream the shock. We would like to note, however, that the 
respective influence of these boundaries discussed in the previous section can 
change if the downstream particle diffusion coefficient is much smaller than the 
upstream one. In general, the boundary allowing for quick particle escape will 
always preferentially act to increase the particle spectrum inclination, while the 
other one will limit in a higher degree the acceleration time scale. Let us also 
note the fact that a downstream boundary further away than approximately 
two diffusive length scales (cf. Figures 3) has almost no influence on the parti- 
cle spectral index at the shock. Therefore physical processes and/or conditions 
behind that distance are not expected to modify the particle distribution at the 
shock. 

The realistic conditions near astrophysical shocks are expected to involve 
the diffusion coefficients depending on particle momentum. The same will hold 
for the particle escape probability, defined by the boundary distance in our 
model. Therefore, the present results can be used only to a rather general eval- 
uation of the acceleration conditions and compatibility of the observed (possibly 
non-power-law) spectra and time scales with the shock dynamics. A detailed 
modelling requires numerical methods and not-frcqucntly available information 
about the local physical conditions and the boundary conditions. 

Another fact of interest for modelling the generation of highest energy par- 
ticles near the spectrum cut-off should be mentioned in this place. The first 
such particles to appear at the shock are those which have not spent too much 
time diffusing far from the shock, in analogy to the case with escape bound- 
aries. Therefore, at a given energy, the time for these particles to appear can 
be substantially shorter than the respective acceleration time scale given in 
Equation (1.2) or (2.24), the one valid for the steady state particle spectrum. 
This fact is visible in analytic time-dependent solutions (e.g. Drury 1991; our 
solution in Appendix A). 
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was partly done during the visit of MO to Max-Planck- Institut fur Radioas- 
tronomie in Bonn. He is grateful to Prof. R. Wiclcbinski and other colleagues 
from the Institute for their hospitality and valuable discussions. The work of 
MO was supported from the KBN grant PB 1117/P3/94/06. RS acknowledges 
partial support by the Deutsche Forschungsgemeinschaft (Schl 201/8-2). 
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Appendix A: Inversion of Equation (2.17) for L\ = 

L 2 = oo 

In the case of no escape boundaries Li = L 2 = oo we obtain from Equa- 
tions (2.13-16) 

g\{s,x,p) = G 1 (s)G 2 (s) , g 2 {s,x,p) = G 3 (s)G 4 {s) , (Al) 

with 

n ( \ 3( 5o f 2x \ ( p x " " 

Gi(S) = —— rrexp' 



where 



I7 2 (i2- 1) 1 \ Pq 
g( P - P o) CXp(AlVT ^ , (A2) 



G 2 (s) = cxp(-A 2 Vl + t 2 s) , (A3) 



r< I \ 3Qo ( 2x \ ( P 

G3(s) = ^M^I) exp feJ U 
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G 4 (s) = cxp(-At\/l + ns) , (-45) 
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The upstream and downstream distribution functions then are 

rt 



fi(t,x,p)= [ duF 2 {u)F 1 (t - u) (A9a) 
Jo 



and 
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f2(t,X,p) 



duF 4 (u)F 3 (t-u) 



(A9b) 



with the respective Laplace inverse transforms 



F 1 (t)=L- 1 (G 1 
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ln(p/p ) + 
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+ exp kk) crfc ( [w=T) Hp/po) + ^u]^ Jl+ ^) } ' 



F 4 (t) = L~ 1 (g 4 ) = 



3^R\n(p/p ) 
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In deriving the four Laplace inversions (A10) we use the tables of Oberhettingcr 
and Badii (1973). Equations (A9) and (A10) represent the full solutions. 
At the shock (x = 0) we obtain 



t 9rii? 2 (ln(p/p )) S 



(AlOc) 



(AlOd) 
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fo(p,t) = .hA x = o,p,t) 



QQoV^IHp/po) ( p\ 



H(p — po) / duu 3,72 exp 
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4( /_ 1} Hp/po) y/n/(t - u) + V(t - + 
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After infinite time Equation (All) approaches the steady-state solution 



(All) 



fo{p,t -> oo) 



9Qo^ln(p/p ) ( p\ 2 <«- 1 > 



f°° ( u 
H(p — p ) / duu~ 3 ^ 2 exp [ 

Jo V T 2 
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3Q H(p - po) ( p \ " 



C/ 2 (fl-l) Vpo 



(A12) 



where we used the integral 



J dxx- 3 / 2 exp (-x - = 2^ 2 K 1/2 (z) Z - 1 / 2 = 2^e~ z /z . (.413) 
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